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QQ ! Abstract 

■ Improving our understanding of the initial conditions and earliest stages of star 
^ formation is crucial to gain insight into the origin of stellar masses, multiple sys- 
Q_i: tems, and protoplanetary disks. We review the properties of low-mass dense cores 
Q as derived from recent millimeter /submillimeter observations of nearby molecu- 
lar clouds and discuss them in the context of various contemporary scenarios for 
^ I cloud core formation and evolution. None of the extreme scenarios can explain all 
I— I, observations. Pure laminar ambipolar diffusion has relatively long growth times 
I ' for typical ionization levels and has difficulty satisfying core lifetime constraints. 
I> ' Purely hydrodynamic pictures have trouble accounting for the inefficiency of core 
. formation and the detailed velocity structure of individual cores. A possible favor- 
', able scenario is a mixed model involving gravitational fragmentation of turbulent 
[ molecular clouds close to magnetic criticality. The evolution of the magnetic field 
and angular momentum in individual cloud cores after the onset of gravitational 
collapse is also discussed. In particular, we stress the importance of radiation- 
, magnetohydrodynamical processes and resistive MHD effects during the proto- 
stellar phase. We also emphasize the role of the formation of the short-lived 
first (protostellar) core in providing a chance for sub-fragmentation into binary 
systems and triggering MHD outflows. Future submillimeter facilities such as 
5^ ', Herschel and ALMA will soon provide major new observational constraints in 
this field. On the theoretical side, an important challenge for the future will be to 
link the formation of molecular clouds and prestellar cores in a coherent picture. 
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1. Introduction: Dense cores and the origin of the IMF 

Stars form from the gravitational collapse of dense cloud cores in the molecular 
interstellar medium of galaxies. Studying and characterizing the properties of 
dense cores is thus of great interest to gain insight into the initial conditions and 
initial stages of the star formation process. 

Our observational understanding of low-mass dense cores has made signifi- 
cant progress in recent years and three broad categories of cores can now be 
distinguished within nearby molecular clouds, which possibly represent an evolu- 
tionary sequence: starless cores, prestellar cores, and "Class 0" protostellar cores. 
Starless cores are possibly transient concentrations of molecular gas and dust 
without embedded young stellar objects (YSOs), typically observed in tracers 
such as C^^O (e.g. Onishi et al. 1998), NH3 (e.g. Jijina, Myers, & Adams 1999), 
or dust extinction (e.g. Alves et al. 2007), and which do not show evidence of 
infall. Prestellar cores are also starless (M^ = 0) but represent a somewhat denser 
and more centrally-concentrated population of cores which are self-gravitating, 
hence unlikely to be transient. They are typically detected in (sub) millimeter 
dust continuum emission and dense molecular gas tracers such as NH3 or N2H+ 
(e.g. Ward-Thompson et al. 1994; Benson & Myers 1989; CaselH et al. 2002), 
often seen in absorption at mid- to far-infrared wavelengths (e.g. Bacmann et al. 
2000, Alves et al. 2001), and frequently exhibit evidence of infall motions (e.g. 
Gregersen & Evans 2000). Conceptually, all prestellar cores are starless but only 
a subset of the starless cores evolve into prestellar cores; the rest are presumably 
"failed" cores that eventually disperse and never form stars. In practice, prestel- 
lar cores are characterized by large density contrasts over the local background 
medium. Specifically, the mean densities of observed prestellar cores exceed the 
mean densities of their parent clouds by a factor ^ 5-10, while their mean col- 
umn densities exceed the background column densities by a factor ^ 2. For 
comparison, a critical self-gravitating Bonnor-Ebert isothermal spheroid has a 
mean density contrast pbe/Pcxi ~ 2.4 (e.g. Lombardi & Bertin 2001) and a mean 
column density contrast SBE/^^cxt ~ 1-5 over the external medium. 

Finally, Class cores/objects are young accreting protostars observed early 
after point mass formation while most of the mass of the system is still in the 
form of a dense core/envelope as opposed to a YSO (M^ <C Menv) (Andre, Ward- 
Thompson, Barsony 1993). They are believed to result from the gravitational col- 
lapse of prestellar cores. Class protostars themselves evolve into Class I objects 
with Menv < -Mt (Lada 1987; Andre & Montmerle 1994) as the protostellar enve- 
lope dissipates through accretion and ejection of circumstellar material. Class I 
objects subsequently evolve into (Class II and Class III) pre-main sequence stars 
surrounded by a circumstellar disk (optically thick and optically thin in the near- 
/mid-IR, respectively), but lacking a dense circumstellar envelope (Menv ~ 0). 

Improving our understanding of the formation and evolution of dense cores in 
molecular clouds is crucially important since there is now good evidence that these 
early stages largely control the origin of the stellar initial mass function (IMF). 
Indeed, observations indicate that the prestellar core mass function resembles the 
IMF (see §[2] below), suggesting that the effective reservoirs of mass required for 
the formation of individual stars are already selected at the prestellar core stage. 
Furthermore, it is at the very end of the prestellar stage that multiple systems are 
believed to form and during the protostellar stage that a fraction of the prestellar 
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Figure 1. SCUBA 850 dust continuum map of the NGC 2068 protocluster extracted 
from the mosaic of Orion B by Motte et al. (2001). A total of 30 compact prestellar 
condensations (marked by crosses), with masses between 0.4 Mq and ~ 4.5 Mq, are 
detected in this 1 pc x 0.7 pc field. 

core mass reservoir is accreted by the central protostellar components as a result 
of the accretion/ejection process. 

2. Link between the prestellar core mass function and the IMF 

Wide-field (sub)mm dust continuum mapping is a powerful tool to take a cen- 
sus of prestellar dense cores and young protostars within star- forming clouds. The 
advent of large-format bolometer arrays on (sub)millimeter radiotelescopes such 
as the IRAM 30m and the JCMT has led to the identification of numerous cold, 
compact condensations that do not obey the Larson (1981) self-similar scaling 
laws of molecular clouds and are intermediate in their properties between diffuse 
CO clumps and infrared young stellar objects (cf. Andre et al. 2000 and Ward- 
Thompson et al. 2007 for reviews). As an example, Fig.[T]shows the condensations 
found by Motte et al. (2001) at 850 //m in the NGC 2068 protocluster (Orion B). 
Such highly concentrated (sub)millimeter continuum condensations are at least 3 
to 6 orders of magnitude denser than typical CO clumps (e.g. Kramer et al. 1998) 
and feature large 50%) mean column density contrasts over their parent back- 
ground clouds, strongly suggesting they are self- gravitating. The latter is directly 
confirmed by line observations in a number of cases. When available, the virial 
masses of the condensations indeed agree within a factor of ~ 2 with the masses 
derived from the (sub)millimeter dust continuum (e.g. Andre et al. 2007). A small 
fraction of these condensations lie at the base of powerful jet-like outflows and 
correspond to Class objects. However, the majority of them are starless/jetless 
and appear to be the immediate prestellar progenitors of individual protostars or 
protostellar systems. 

In particular, as first pointed out by Motte, Andre, Neri (1998) in the case 
of the p Ophiuchi (L1688) cloud, the mass distribution of these starless dust 
continuum condensations is remarkably similar in shape to the stellar IMF (see 
Fig. ED- This was consistently found by a number of independent groups in the 
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Figure 2. Differential (dN/dlogM) mass distribution of the 229 starless dust continuum 
condensations detected at 850 /zm with SCUBA in the Orion A/B cloud complex exclud- 
ing the crowded OMCl and NGC 2024 regions (histogram with error bars - from Motte 
et al. 2001, Johnstone et al. 2001, and Nutter & Ward-Thompson 2007). This prestellar 
core sample is estimated to be complete down to ~ 0.3 Mq. A two-segment power law fit 
and a lognormal fit are shown for comparison. The lognormal fit peaks at ~ 1.1 Mq and 
has a standard deviation of ^ 0.41 in logioM. For reference, the dash-dotted curve shows 
the shape of the single-star IMF (e.g. Kroupa 2001) and the dashed curve corresponds 
to the IMF of multiple systems (e.g. Chabrier 2005). (The lognormal part of the latter 
peaks at 0.25 Mq and has a standard deviation of 0.55 in logioM.) The dotted line shows 
a dN/dlogM (x M~°-^ power-law distribution corresponding to the typical mass spectrum 
found for low-density CO clumps (see Blitz 1993 and Kramer et al. 1998). 

past few years (e.g. Testi &: Sargent 1998; Johnstone et al. 2000, 2001; Motte 
et al. 2001; Stanke et al. 2006; Enoch et al. 2006; Nutter & Ward-Thompson 
2007) in nearby star-forming regions such as p Ophiuchi, Serpens, Orion A & 
B, and Perseus. In all of these clouds, the observed prestellar core mass function 
(CMF) is consistent with the Salpeter (1955) power-law IMF at the high-mass end 
(dN/dlogM oc M~^'^^), and thus significantly steeper than the mass distribution 
of diffuse CO clumps (dN/dlogM oc M^*^'^ over at least three decades in mass 
- e.g. Blitz 1993; Kramer et al. 1998). The difference presumably arises because 
CO clumps are primarily structured by supersonic turbulence (e.g. Elmegreen 
& Falgarone 1996) while prestellar condensations are largely free of supersonic 
turbulence and clearly shaped by self-gravity (e.g. Motte et al. 2001, Andre et 
al. 2007 and § 13.21 below). The slope of the observed CMF becomes shallower 
than the Salpeter power law and more similar to the slope of the typical CO 
clump mass distribution at the low-mass end. Based on the results of present 
core surveys, the entire prestellar CMF can generally be fit equally well with 
either a two-segment broken power law or a lognormal distribution down to the 
completeness limit of the observations. This is illustrated in Fig. [2] which shows 
the differential CMF in log-log format for a sample of 229 starless submillimeter 
continuum cores in the Orion A/B complex identified in the SCUBA surveys of 
Motte et al. (2001) and Johnstone et al. (2001, 2006). 

Note that there is some discussion in the literature (e.g. Reid &: Wilson 2006) 
as to whether the differential or the cumulative form of the CMF should be used. 
The differential form is more intuitive to interpret as it can be fit using the 
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standard least-squares technique after assigning simple Poisson error bars. It is 
however inadequate when the total number of objects is small ( ^ 100) as it is 
significantly affected by the arbitrary choice of mass bins. The cumulative form, 
which is independent of binning, is preferable when dealing with small samples, 
but some care must be taken when comparing it with model distributions based on 
least-squares fits (cf. Reid & Wilson 2006). In practice, using the non-parametric 
Kolmogorov-Smirnov (K-S) test to compare the cumulative form of the CMF 
with model distributions is more reliable than least-squares fitting. The K-S test 
confirms that, within statistical uncertainties, the prestellar CMF observed in 
the above-mentioned nearby clouds is indistinguishable in shape from the stellar 
IMF, although the significance remains limited by the relatively small size of 
present prestellar core samples. The K-S test also shows that the observed CMF 
differs from the shallow power-law mass distribution of CO clumps at a very 
high significance level (for instance, in the Orion case illustrated in Fig. [2l the 
probability that the two mass distributions are statistically similar in shape is 
only P w 2 X lO'^). 

In addition, the median prestellar core mass observed in regions such as p Ophi- 
uchi and Orion (~ 0.2 — 1.5 Mq) is only slightly larger than the characteristic 
~ 0.5 Mq set by the peak of the IMF in dN/dlogM format. Such a close resem- 
blance of the CMF to the IMF in both shape and mass scale is consistent with 
the view that the prestellar condensations identified in (sub)millimeter dust con- 
tinuum surveys are about to form stars on a one-to-one basis, with a fixed and 
relatively high local efficiency, i.e., ecorc = Mi,/Mcore ~ 30 — 50%. 

Interestingly, in a recent near-IR extinction imaging study of the Pipe dark 
cloud, Alves et al. (2007) found a population of 159 starless cores whose mass 
distribution similarly follows the shape of the IMF. This finding is reminiscent 
of the CMF results obtained from the (sub) millimeter dust continuum, although 
it is important to stress that most of the starless cores in the Pipe Nebula are 
gravitationally unbound objects confined by external pressure (Lada et al. 2008). 
Hence, they do not qualify as prestellar cores and a large fraction of them may 
never evolve into stars. Assuming nevertheless that most of them will evolve into 
self-gravitating prestellar cores and subsequently collapse into stars, the Alves et 
al. (2007) result suggests that the IMF may be determined even earlier than the 
prestellar stage. 

Appealing as a direct connection between the prestellar CMF and the IMF 
might be, several caveats should be kept in mind. First, although core mass 
estimates based on optically thin (sub)millimeter dust continuum emission are 
straightforward, they rely on uncertain assumptions about the dust (tempera- 
ture and emissivity) properties (e.g. Stamatellos et al. 2007a). Second, current 
determinations of the CMF are limited by small-number statistics in any given 
cloud and may be affected by incompleteness at the low-mass end (e.g. John- 
stone et al. 2000). With Herschel, the future submillimeter space telescope to 
be launched in 2008, it will be possible to dramatically improve on the statistics 
and to largely eliminate the mass uncertainties through direct measurements of 
the dust temperatures (cf. Andre &: Saraceno 2005). Third, in some regions such 
as the p Ophiuchi cloud, the shape of the CMF is in better agreement with the 
IMF of individual field stars than with the IMF of multiple systems (e.g. Andre 
et al. 2007). This is surprising since the spatial resolution of current surveys for 
prestellar cores (~ 2000 AU at best) is not sufficient to probe core multiplicity. 
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Furthermore, multiple systems are believed to form after the prestellar stage by 
subsequent dynamical fragmentation during the collapse phase, close to the time 
of protostar formation (e.g. Goodwin et al. 2007). Thus, one would expect the 
masses of prestellar cores to be more directly related to the masses of multi- 
ple systems than to the masses of individual stars. It is possible that a fraction 
at least of the cores observed at masses lower than the peak of the CMF in 
dN/dlogM format (e.g., ~ 1.1 Mq in Fig. [2|) are not gravitationally bound, hence 
not prestellar in nature (cf. Andre et al. 2007). 

Last but not least, there is a potential timescale problem. As pointed out by 
Clark, Klessen, &: Bonnell (2007), if the lifetime of prestellar cores depends on 
their mass, then the observed mass distribution is not necessarily representative 
of the intrinsic CMF (see also Elmegreen 2000). This is due to the fact that an 
observer is more likely to detect long-lived cores than short-lived cores. In prac- 
tice, however, the mean densities of prestellar cores are essentially uncorrelated 
with their masses, so that there is no systematic dependence of the dynamical 
timescale on the mass. The importance of the potential timescale bias can be as- 
sessed by considering a weighted core mass function in which each core is assigned 
a weight equal to = -^^iTsy (instead of 1), where (tg) is the average free-fall 
time of the sampled cores. Such a weighting makes it possible to recover the 
intrinsic shape of the CMF if the lifetime of each core is proportional to its free- 
fall time. Andre et al. (2007) applied this technique to the sample of 57 starless 
dust continuum condensations identified by Motte et al. (1998) at 1.2 mm in the 
p Ophiuchi cloud. In this case, the above-mentioned weighting does not change 
the high-mass end of the CMF and only affects the low-mass end: it renders the 
entire p Oph CMF derived above the ~ 0.1 completeness level remarkably 
consistent with a single, Salpeter power-law mass distribution (see Fig. 8 of Andre 
et al. 2007). We conclude that the steep, Salpeter-like slope of the CMF at the 
high-mass end is robust, but that the departure from a single power-law distri- 
bution at the low-mass end, in the form of a break near the median core mass 
(cf. Fig. [2]), is less robust. 

Despite these limitations, the observational findings summarized in this section 
are very encouraging as they support scenarios according to which the bulk of the 
IMF is partly determined by pre-collapse cloud fragmentation (e.g. Larson 1985, 
2005; Elmegreen 1997; Padoan & Nordlund 2002). The finding that the high-mass 
end of the prestellar CMF is substantially steeper than the dN/dlogM oc M"*^'^ 
mass distribution of low-density CO clumps is very significant: while most of the 
mass of the low-density CO medium is contained in the largest, most massive 
clumps, most of the prestellar mass destined to evolve into stars is in small, low- 
mass cores. Clearly, one of the keys to the problem of the origin of the IMF lies 
in a good understanding of the processes responsible for the formation of prestel- 
lar cores/condensations out of low-density structures within molecular clouds. 
However, it is likely that additional processes, such as subfragmentation into bi- 
nary/multiple systems and more generally mechanisms controlling the value of 
the star formation efficiency at the core level (ecore)) also play an important role 
and, in particular, are required to generate the low-mass (M < 0.3 Mq) end of 
the IMF (cf. Bate et al. 2003, Ballesteros-Paredes et al. 2006). In the following 
sections, we discuss core formation and core subfragmentation models in turn. 
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3. Core formation models vs. observational constraints 

The mechanisms by which presteUar cores form and evolve in molecular clouds 
are the subject of a major theoretical debate at the moment. There is little doubt 
that self-gravity ultimately plays a dominant role and it has even been proposed 
that dense cores may form by purely gravitational fragmentation (e.g. Larson 
1985; Hartmann 2002). However, the respective roles of magnetic fields and in- 
terstellar turbulence in regulating the core/star formation process are highly con- 
troversial. In particular, the classical picture of slow, quasi-static core formation 
by ambipolar diffusion in magnetically-supported clouds (e.g. Mouschovias 1987; 
Shu et al. 1987, 2004) has been seriously challenged by a new, much more dy- 
namic paradigm, which emphasizes the role of supersonic turbulence in support- 
ing clouds on large scales and generating density fluctuations on small scales 
(e.g. Padoan & Nordlund 2002; Mac Low & Klessen 2004). Conceptually, core 
formation models may be conveniently divided up into four categories, depend- 
ing on whether linear or non-linear (turbulent) perturbations initiate core for- 
mation, and on whether magnetic fields are dynamically dominant or not (see 
Table [1] below) . In this classification, the "standard" picture and the new tur- 
bulent paradigm correspond to two extreme models, according to which cores 
form by magnetically-regulated gravitational fragmentation and super-Alfvenic 
turbulent fragmentation, respectively. In all turbulent models, cores are initially 
formed by cloud material compressed by shocks arising from supersonic turbu- 
lence. There are however several versions of the dynamic picture of core/star 
formation which mainly differ in the way they explain the origin of the IMF. In 
the Padoan &; Nordlund (2002) scenario, the IMF is almost entirely set by the 
properties of interstellar turbulence at the prestellar stage, while in the alterna- 
tive model proposed by Bate & Bonnell (2005), turbulence is largely irrelevant for 
the IMF which originates from competitive accretion and dynamical interactions 
at the protostellar stage. The Klessen &: Burkert (2000) scenario is intermediate 
between these two extremes in that both turbulence and dynamical interactions 
play a role in shaping the IMF. The various turbulent models also differ depend- 
ing on whether the turbulence is freely decaying (e.g. Tilley & Pudritz 2007) or 
continuously driven (e.g. Vazquez-Semadeni et al. 2005). 



3.1. Theoretical description of cloud fragmentation models 

The classical problems of fragmentation of a sheet-like layer or a cylinder are 
relevant to star formation because they show that there is a preferred scale of 
gravitational fragmentation as soon as one considers a structured (i.e., flattened 
or filamentary) region with some scale length H (see Larson 1985). This is not 
the case for a uniform medium, in which case the Jeans analysis shows that 
the largest possible scale has the fastest growth rate. Sheets and filaments can 
be easily generated by the formation process of the molecular cloud itself or by 
subsequent internal turbulent or gravitational motions. For isothermal sheets, 
H = Cg/(7rGS), where Cg is the isothermal sound speed, and S is the column 
density of the sheet. For highly flattened sheets, the preferred fragmentation 
scale is Xm = ^irH, while for a layer with the extended isothermal atmosphere 
calculated by Spitzer (1942), Am = ^A'kH (Simon 1965). These two length scales 
likely bracket the possibilities for more realistic isothermal non-magnetic sheet- 
like configurations. The growth time of the fragmentation instability is essentially 
the dynamical timescale t^ ~ H/cs for the cases described above, and is more 
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generally identified with the free-fall timescale ~ Xj \JGp that applies to both 
unpressured sheets and those confined by a strong external pressure. 

However, purely gravitational fragmentation instability of an entire molecular 
cloud on the dynamical timescale given by the mean column density and tempera- 
ture is ruled out (Zuckerman and Palmer 1974), due to the observed low efficiency 
of Galactic star formation. Nevertheless, it remains a relevant concept to under- 
standing fragmentation in the cluster- forming subregions of molecular clouds. 
A proposed explanation for the low efficiency of star formation is the so-called 
"standard model" of star formation (Shu, Adams, & Lizano 1987; Mouschovias 
1987), in which cores are formed on a diffusive timescale (much longer than than 
a dynamical timescale), due to the ambipolar drift of neutrals past dynamically 
dominant magnetic fields. This mode of star formation is often erroneously labeled 
as "isolated" star formation. In reality, ambipolar diffusion driven core formation 
is also a fragmentation process as surely as its nonmagnetic counterpart - the 
difference is primarily in the timescale of evolution. One may however make the 
following distinction based on the mass-to- magnetic flux ratio, M/<I>, compared to 
the critical value, (M/<I>)^j.j^., necessary for support against gravitational collapse. 
Molecular cloud envelopes may be magnetically dominated, with subcritical to 
transcritical mass-to-flux ratios, i.e., //env = (M/<I>) / (M/$)^j.;^. ^ 1. Cloud en- 
velopes also have low mean density and relatively high levels of ionization, and 
thereby evolve on such a long timescale for fragmentation that no runaway has 
essentially occurred by the time we observe the clouds. Conversely, the cluster- 
forming cores (and also regions of weak clustering found in the Taurus molecular 
cloud) may be magnetically supercritical with //dus > 1 and also have greater 
mean density, resulting in fragment formation on relatively shorter timescales 
and lengthscales. 

For the above reasons and because Zeeman measurements (see the compilation 
of Crutcher 1999 and ^ l3.2.3l below) establish that mass-to-flux ratios are clustered 
about the critical value, it is important to study cloud fragmentation including 
the effect of dynamically significant magnetic fields and also ambipolar diffusion, 
if possible. Magnetized sheets can be studied most easily in two limits. If the 
formation process of the sheet is by dynamical compression, say from stellar 
winds or supernovae, the ambient interstellar magnetic field may be swept up 
into the expanding shell and be oriented primarily in the plane of the sheet. 
Conversely, if self-gravity plays an important role in the formation of the sheet, 
then it will be flattened along the mean direction of the magnetic field. 

When the magnetic field is in the plane of the sheet, the character of fragmenta- 
tion and the fate of the cloud can be further categorized into two subcases (Nagai, 
Inutsuka, & Miyama 1998). If the geometrical thickness of the sheet is comparable 
to (or larger than) the natural scale height (~ Cg/ y/Gpc)-, or equivalently, if the 
ambient pressure is much smaller than the midplane pressure, compressional mo- 
tions along the magnetic field lines result in the formation of filamentary clouds 
elongated perpendicular to the field lines. The line-mass (mass per unit length) of 
the filament is larger than that of the (isothermal) equilibrium filament so that 
collapse toward the axis of the filament continues until temperature increases 
(Inutsuka & Miyama 1992). In this case, the characteristic (i.e., minimum) mass 
scale for fragmentation is determined by the final fragmentation of the filamentary 
cloud, and significantly smaller than the initial Jeans mass of the sheet (Inutsuka 
& Miyama 1997). The resulting mass scale can be called "the minimum Jeans 
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Table 1. Summary of main features of core formation models. 



Core Formation Scenarios 



Models 



Gravitational Fragmentation 
(linear perturbations) 



Turbulent Fragmentation 
(non-linear perturbations) 



Weak B 



Short (few Myr) timescale 
infall mildly supersonic 
ordered curved field lines 
initial CMF very narrow 



Very short (< Myr) timescale 
infall highly supersonic 
field lines distorted 
initial CMF is broad, IMF-like 



Strong B 



Long (~ 10 Myr) timescale 
subsonic infall 

small field line curvature 
initial CMF very narrow 



Short (few Myr) timescale 
subsonic relative infall and 
supersonic systematic speeds 
ordered field lines 
initial CMF is broad, IMF-like 



mass" and its dependence on the initial temperature and metallicity of the cloud 
was discussed by Masunaga & Inutsuka (1999). The expected core mass function 
was analytically derived by Inutsuka (2001) using the Press-Schechter formalism. 
The overall timescale of the whole process is on the order of the initial free-fall 
time (% ~ l/^/Gp'c)- 

In contrast, if the sheet-like cloud is confined by strong external pressure due to 
an ambient warm (or hot) medium, the thickness of the sheet is smaller than 
~ c^l \/Gpc- The timescale for fragmentation is still on the order of the gravita- 
tional free-fall time iff, but sound waves can now propagate many times across 
the (thin) sheet within this timescale. Thus, the fragmentation is in an incom- 
pressible mode and the axes of the resulting filaments are parallel to the direction 
of the magnetic field. In this case, the filaments may fragment into cores whose 
masses can be smaller than the Jeans mass. In effect, this second case provides a 
mechanism for generating gravitationally stable cores (Inutsuka & Miyama 1997). 

When the magnetic field is oriented perpendicular to the sheet, recent sim- 
ulations of cloud fragmentation which include ambipolar diffusion have allowed 
an extensive parameter study of the interplay of magnetic fields, ambipolar dif- 
fusion, and turbulence in the formation of cores (Basu & Ciolek 2004; Basu, 
Ciolek, &; Wurster 2008). Given the standard cosmic-ray induced ionization frac- 
tion Xi ~ 10^'^(n/10^ cm^^)^"*^/^ (Elmegreen 1979; Nakano 1979), there is an 
observationally distinguishable difference between subcritical and supercritical 
fragmentation. Subcritical fragmentation leads to subsonic infall motions onto 
cores and within them, while supercritical cloud fragmentation leads to extended 
supersonic infall on the core scale ~ 0.1 pc and even beyond. 

Figure [3] shows a comparison of column density structure and velocity vectors 
for four different models (from Basu, Ciolek, & Wurster 2008; Basu et al. 2008). 
Each image is shown at the time of runaway collapse of the first supercritical 
dense core that forms in the simulation, with a column density enhancement of a 
factor of 10. The top two images are for subcritical (left) and supercritical (right) 
clouds evolving from linear initial perturbations. The bottom two images are for 
subcritical (left) and supercritical (right) clouds evolving from highly nonlinear 
(turbulent) initial conditions. Each image represents a physically distinct path to 
core/star formation. Velocity vectors are overlaid, although the normalizations 
may differ (see caption). The subcritical models have /io = 0.5 (i.e., the mass- 
to-flux ratio is half the critical value) and the supercritical models have ^uq = 2. 
Amongst the gravitational fragmentation models (linear initial perturbations). 
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Figure 3. Image and contours of column density and velocity vectors of neutrals, for 
four different models at the time of runaway collapse of the first supercritical dense core 
(Basu, Ciolek, & Wurster 2008; Basu et al. 2008). Top left: gravitational fragmentation 
for a subcritical model /ig = 0.5. Top right: the same for a supercritical model with 
fiQ = 2.0. Bottom left: turbulent fragmentation for = 0-5- Bottom right: turbulent 
fragmentation for = 2.0. The color table is applied to the logarithm of column density 
and the contour lines represent values of column density enhancement in multiplicative 
increments of 2^/^, having the values [0.7,1.0,1.4,2,2.8,4.0,...]. The horizontal or vertical 
distance between tips of velocity vectors corresponds to a speed 0.5 Cg in the top panels, to 
1.0 Cg in the bottom left panel, and to 3.0 Cg in the bottom right panel. Spatial coordinates 
are normalized to 27r7J, the wavelength of maximum growth rate in the limit of no 
magnetic field and external pressure. 



there is extended and mildly supersonic motions only in the supercritical model, 
vi^hereas the subcritical model has a maximum speed of only ~ 0.4cs. For the 
subcritical gravitational fragmentation model the runaway collapse to the first 
core occurs at a time 7.4 Myr, whereas for the supercritical model it occurs at 0.83 
Myr, assuming a background column density 10^^ cm~^ and temperature 10 K. 
Fragmentation of a transcritical (/io — 1) cloud is qualitatively unique in that the 
fragmentation spacing is much larger than the typical value ~ lirH that is valid 
for both highly supercritical and highly subcritical clouds (Ciolek &: Basu 2006). 
The timescale of transcritical fragmentation is intermediate to the subcritical 
and supercritical cases, but closer to the former. Here, we focus primarily on the 
distinction between decidedly subcritical and supercritical clouds. The bottom 
panels of Fig. [3]show the corresponding fiQ = 0.5 (left) and /io = 2 (right) models, 
now with turbulent initial perturbations with rms amplitude Va = 2cs in each of 
Vx and Vy. The power spectrum is such that cx k""^, so most of the energy 
is in the largest scale modes: the result is immediate large-scale compressive 
motions in the cloud. The = 0.5 model undergoes rapid ambipolar diffusion 
during the first compression, but still rebounds due to stored magnetic energy. It 
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Figure 4. Histograms of masses contained within regions with column density enhance- 
ment above a factor of two, measured at the end point of simulations with /io = 2.0. The 
left panel corresponds to gravitational fragmentation and the right panel corresponds to 
turbulent fragmentation. See text for details. Each figure is obtained from the compilation 
of results of a large number of simulations. The bin width is 0.1. 

subsequently oscillates several times before continuing ambipolar diffusion causes 
runaway collapse in the highest density region. The total time for this process is 
1.2 Myr, similar to that of supercritical gravitational fragmentation. In contrast, 
the /^o = 2 model goes into prompt collapse during the first compression, on 
a timescale of merely 4.0 x 10'^ yr = 0.04 Myr. Systematic supersonic motions 
exist in both turbulent models, but the relative infall is generally subsonic in the 
subcritical model. A comparison of the two bottom panel images shows that the 
supercritical model is extremely filamentary, while the subcritical model is much 
less so since it has had a chance to rebound from the initial compression. These 
turbulent simulations are done in the thin-disk approximation and are consistent 
with the earlier results of Li &: Nakamura (2004) and Nakamura & Li (2005), also 
done using the same approximation. Recent fully three-dimensional simulations 
(Kudoh et al. 2007; Kudoh & Basu 2008) also confirm the results. A summary 
of main outcomes of the four different modes of core/star formation is given in 
Tabled! 

The results described above are used to generate statistics of core masses, 
sizes, shapes, etc. Each model with a unique set of parameters is run ~ 100 times 
in order to generate a meaningful measure of the various outcomes arising from 
different realizations of the initial random perturbations. Thresholding techniques 
are used to identify cores - for details of the technique, see Basu, Ciolek, & 
Wurster (2008). Figure [3] shows a comparison of the initial core mass function 
(CMF) for supercritical models with /xq = 2 but without and with turbulent 
initial perturbations, respectively. We emphasize these are initial CMF's because 
they refiect the status of cores at the time that the first supercritical dense core 
undergoes runaway collapse. The clear distinction to be made between the two 
models is that the linear initial perturbation case (left panel) has an extremely 
sharp peak (consistent with a preferred mass scale for fragmentation in the linear 
theory), although there is a broader tail at the low mass end due to some cores 
that are very young and just emerging above the threshold. In contrast, the 
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Figure 5. SCUBA 850 fim dust continuum map (colorscale) of the Western part of the 
Perseus cloud complex. Contours of C^'*O(l-0) integrated intensity are overlaid. The black 
boundary indicates the total area mapped with SCUBA. Note that dust continuum cores 
are detected only localized sub-regions, such as NGC 1333, within the complex. (From 
Hatchell et al. 2005) 

turbulent initial condition case develops a broad tail of high mass cores, due to 
much of the turbulent power being on large scales. However, the high mass cores 
often have quite disturbed structures (see Fig. [3] lower panels) and it is not clear 
that they would subsequently collapse monolithically. 

The narrow initial CMF seems to be a problem for purely gravitational frag- 
mentation models. However, a broader distribution may yet be possible if the 
newly formed cores continue to accrete from their environment. If the accretion 
from the environment is nonuniform from core to core, an initially narrow and/or 
lognormal CMF may become skewed so that a high mass tail develops. For exam- 
ple, Basu & Jones (2004) have generalized an earlier result of Myers (2000) and 
shown that an initially lognormal CMF and an exponential distribution of subse- 
quent accretion times from the cloud results in a CMF that is like a lognormal at 
low masses but has a power-law tail at high masses. The model formally requires 
that the accretion rate onto a core is linearly proportional to the instantaneous 
core mass. Additionally, Basu & Jones (2004) show that a different accretion law 
can also lead to a broadened near-power-law tail. 

3.2. Observational diagnostics 

In principle, it should be possible to discriminate between these various core 
formation scenarios based on detailed observational studies of the characteristics 
of prestellar cores and young protostars. 

3.2.1. Core formation efficiency and spatial distribution from surveys 

A number of very large submillimeter continuum surveys of nearby cloud com- 
plexes have been completed recently which provide the spatial distribution of 
cores within these complexes and set constraints on the efficiency of the core 
formation process (e.g. Johnstone et al. 2004, Hatchell et al. 2005, Enoch et al. 
2006, Motte et al. 2007). As an illustration. Figure [5] shows part of the 3 deg^ 
SCUBA 850 ^m survey of the Perseus cloud complex by Hatchell et al. (2005). 
Such extensive surveys show that prestellar cores and Class protostars are found 
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Figure 6. Plot of inferred core lifetime against mean volume density for various samples 
of starless and prestellar cores. The two dashed lines correspond to one and ten free-fall 
times, respectively. (From Jessop & Ward- Thompson 2000 and Kirk et al. 2005.) 

in localized sub-regions within molecular clouds which occupy only a very small 
fraction of their volume. These localized active sub-regions often correspond to 
cluster-forming clumps associated with embedded near-IR clusters. On this ba- 
sis, it has been suggested that there may be a threshold in background column 
density (or equivalently visual extinction at Ay ~ 5 — 10) for core formation 
(Onishi et al. 1998, Johnstone et al. 2004). Observationally, however, establish- 
ing the presence (or absence) of such a threshold is difficult since there are also 
detection thresholds. According to Hatchell et al. (2005), there is no real thresh- 
old but the probability of forming a prestellar core is a steeply rising function 
of background column density. In any case, the results of these wide-field sur- 
veys for cores clearly demonstrate the global inefficiency of the core formation 
process: the fraction of cloud mass observed in the form of prestellar cores is 
very low ( < 1 - 20% - Hatchell et al. 2005, Nutter et al. 2006). Furthermore, 
there is some evidence that core/star formation in the observed cluster-forming 
clumps has been induced by external triggers (Nutter et al. 2006, Kirk et al. 
2006). These results are broadly consistent with the view that the low-density 
envelopes of molecular clouds are supported against collapse by magnetic fields, 
as in the classical ambipolar diffusion picture (Mouschovias 1987, Shu et al. 1987, 
McKee 1989). 

3.2.2. Core lifetimes 

In the turbulent paradigm, cloud cores are always dynamically evolving (e.g. 
Dib et al. 2007) and survive for at most a few free-fall times. 
Observationally, a rough estimate of the lifetime of starless cores can be obtained 
from the number ratio of cores with and without embedded YSOs in a given core 
sample. Using this technique, Lee &: Myers (1999) found that the typical lifetime 
of starless cores with average volume density ~ 10^ cm~^ was ~ 1 — 1.5 x 10^ yr. 
Furthermore, by considering several samples of isolated cores spanning a range 
of core densities, Jessop & Ward-Thompson (2000) established that the typical 
core lifetime decreased as the mean volume density in the core sample increased 
(Fig. E]- see also Kirk et al. 2005). As can be seen in Fig. El all of the observed 
lifetimes lie between one free-fall time, which is the timescale expected in free- 
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Figure 7. Plot of observed mass-to-magnetic flux ratio, in units of the critical value and 
divided by 3 to correct for projection bias, against cloud/core column density. Dots are 
for Zeeman detections of By (above 3cr); stars are for Chandraskehar- Fermi estimates of 
Bj^; open triangles are lower limits (corresponding to observed upper limits to the field 
strength). (From Heiles & Crutcher 2005 - see also Crutcher 2004 for details.) 

fall collapse, and ten free-fall times, the timescale expected for highly subcritical 
cores undergoing ambipolar diffusion. The observed timescales are typically longer 
than the free-fall time by a factor of ~2-5 in the density range of 10^-10^ cm^^. 
This suggests that starless cores cannot all be rapidly evolving, at variance with 
purely hydrodynamic scenarios of core formation (e.g. Ballesteros-Paredes et al. 
2003), but in agreement with numerical simulations of moderately supercritical, 
turbulent molecular clouds (Galvan-Madrid et al. 2007). The observed timescales 
also appear to be too short for all of the cores to form from a highly subcritical 
state. These statistical estimates of core lifetimes are however quite uncertain 
since they assume that all of the observed cores follow the same evolutionary 
path and that the core/star formation rate is constant. 

3.2.3. Magnetic field measurements 

In principle, observations of magnetic fields can provide a strong discriminator 
between the two extreme paradigms of the core/star formation process. Two main 
observational techniques have been used to estimate the magnetic field strength 
in cloud cores. First, the Zeeman effect in, e.g., the 18-cm lines of OH, gives 
a direct measurement of the line-of-sight component of the magnetic field, but 
relatively few positive detections have been obtained. For existing detections, 
a very good correlation is found between the magnetic field strength Bio^ and 
n^'^ ay, where n is the gas density and is the line-of-sight velocity dispersion 
observed in optically thin molecular line tracers (Crutcher 1999, Basu 2000). 
This shows that magnetic fields play a dynamically important role during the 
contraction of at least some cloud cores and is consistent with cloud turbulence 
being MHD in character. Second, maps of linearly polarized dust emission at 
submillimeter wavelengths, combined with the Chandrasekhar and Fermi method, 
provide an indirect estimate of the plane-of-the-sky component of the magnetic 
field (e.g. Ward-Thompson et al. 2000, Crutcher et al. 2004). The available B- 
field measurements based on either the Zeeman technique or the Chandrasekhar 
and Fermi method are summarized in Fig. [7] which plots the inferred mass-to- 
magnetic flux ratio corrected for projection effects against the column density of 
each core (see Crutcher 2004 and Heiles &: Crutcher 2005). It can be seen that 
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all cloud cores are scattered around the critical mass-to-flux ratio in this plot, 
suggesting that, on average, cores are close to magnetically critical. Furthermore, 
there is some hint in Fig. [7] that the mass-to-flux ratio may be systematically 
subcritical at column densities lower than ~ 3 x 10^^ cm~^. This tentative trend 
is weak, however, as it relies only on the locations of Zeeman non-detections in 
the diagram. Thus, existing magnetic fleld measurements do not allow a deflnite 
conclusion to be drawn, even if they seem to favor pictures of the core formation 
process in which the magnetic fleld does play an important role. 

3.2.4. Radial density structure 

The density proflles of isolated prestellar cores are now fairly well known. Two 
methods have been used: (1) mapping the optically thin (sub)millimeter con- 
tinuum emission from the cold dust contained in the cores, and (2) mapping 
the same cold core dust in absorption against the background infrared emission 
(originating from warm cloud dust or remote stars). 

Ward-Thompson et al. (1994, 1999) and Andre et al. (1996) employed the flrst 
approach to probe the structure of prestellar cores (see also Shirley et al. 2000). 
Under the simplifying assumption of spatially uniform dust temperature and 
emissivity properties, they concluded that the radial density proflles of isolated 
prestellar cores were flatter than p(r) oc r^^ in their inner regions (for r ^ Rfiat), 
and approached p(r) oc only beyond a typical radius Rfiat ~ 2500-5000 AU. 

More recently, the use of the absorption approach, both in the mid-IR from 
space (e.g. Bacmann et al. 2000) and in the near-IR from the ground (e.g. Alves 
et al. 2001), made it possible to conflrm and extend the (sub)millimeter emission 
results, essentially independently of any assumption about the dust temperature 
distribution. In some cases, such as L1689B, the absorption studies indicate that 
isolated prestellar cores feature sharp edges defining outer radii Rout ~ 0.1 pc 
(cf. Bacmann et al. 2000). 

The circularly-averaged column density profiles can often be fit remarkably well 
with models of pressure-bounded Bonnor-Ebert spheres, as flrst demonstrated 
by Alves et al. (2001) for B68. This is also the case of cores detected in the 
submillimeter continuum such as L1689B (e.g. Kirk et al. 2005). The quality of 
the flts shows that equilibrium Bonnor-Ebert spheroids provide a good, flrst order 
model for the structure of isolated prestellar cores. In detail, however, there are 
problems with this model. First, the inferred density contrasts (from center to 
edge) are generally larger (i.e., ^ 20-80 - cf. Bacmann et al. 2000 and Kirk et 
al. 2005) than the maximum contrast of ~ 14 for stable Bonnor-Ebert spheres. 
Second, the effective gas temperature needed in the flts is often signiflcantly 
larger than measured core temperatures (e.g. Ward-Thompson et al. 2002, Lai 
et al. 2003). These arguments suggest that prestellar cores are either already 
contracting (see Lee, Myers, Tafalla 2001 and ^ l3.2.5l below) or experiencing extra 
support from static or turbulent magnetic flelds (e.g. Curry & McKee 2000). As 
shown by Bacmann et al. (2000), one way to account for large density contrasts 
and high effective temperatures is to consider models of cores initially supported 
by a static magnetic fleld and evolving through ambipolar diffusion (e.g. Ciolek 
& Mouschovias 1994, Basu & Mouschovias 1994). 

However, good Bonnor-Ebert flts can often be found for dynamically evolv- 
ing "cores" produced by turbulent compression (Ballesteros-Paredes et al. 2003, 
Gomez et al. 2007). Thus, the observed density proflles do not provide a very 
strong diagnostic of proposed core formation models. 
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Figure 8. a) Greyscale image of the line-of-sight velocity dispersion (in units of the 
isothermal sound speed) derived from deep ^^CO(l-O) line observations toward and 
around the prestellar core L1689B (from Andre, Pety, & Bacmann, in prep.). Column 
density contours from Bacmann et al. (2000) are superimposed and delineate the core 
boundaries. The maximum observed value of aios is only ~ 0.34 km/s or 1.7 b) 
Greyscale image of the line-of-sight velocity dispersion superimposed on column den- 
sity contours for a model core obtained in SPH numerical simulations of gravoturbulent 
fragmentation (from Klessen et al. 2005). Note that these simulations produce localized 
maxima where aios 3cs at the shock positions. 

3.2.5. Velocity structure 

Observing the velocity field within and around dense cores is probably the 
most effective way to discriminate between quasi-static and dynamic core forma- 
tion scenarios. If cores are produced by shocks in large-scale supersonic flows, 
large-velocity gradients and local maxima of the line-of-sight velocity dispersion 
are expected in the immediate vicinity of cores (Ballesteros-Paredes et al. 2003; 
Klessen et al. 2005). A much more quiescent ambient velocity field is expected 
in the magnetically-controlled picture (e.g. Nakamura &: Li 2005 - see also § 13.11 
and Fig.E]). 

Briefly, isolated starless cores are characterized by subsonic levels of internal 
turbulence (e.g. Myers 1983, Goodman et al. 1998, Caselli et al. 2002), small ro- 
tational velocity gradients (e.g. Goodman et al. 1993, Caselli et al. 2002), and ex- 
tended, subsonic infall motions (e.g. Tafalla et al. 1998, Lee et al. 2001). Moreover, 
the environment of isolated prestellar cores is also extremely quiescent. This has 
been shown recently through deep mapping of several cores in low-density molec- 
ular gas tracers such as ^^CO(1~0) and C^^O(l-O). Figure [8^ shows a greyscale 
map of the line-of-sight velocity dispersion aios observed in ^'^CO(l-O) toward the 
prestellar core L1689B in the Ophiuchus complex (Andre et al. in prep.). It can 
be seen that aios remains at most transonic {aios < 2 c^, where Cs is the isothermal 
sound speed) everywhere in a region of more than 0.25 pc in diameter around the 
column density peak. Such a quiescent velocity field is at variance with purely 
hydrodynamic models of gravoturbulent fragmentation (e.g. Klessen et al. 2005). 
While these models successfully produce a fair amount (~ 25%) of cores with 
subsonic internal velocity dispersions, corresponding to dense, post-shock stag- 
nation points at the intersection of converging flows, they also produce highly 
supersonic maxima of aios in the low column density gas around the cores (cf. 
Fig. [8)3), which are not observed. Current observations therefore provide strong 
indirect evidence that the evolution of isolated prestellar cores is magnetically 
controlled. 

The environment of individual cores in cluster-forming regions is known to be 
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more turbulent (e.g. Caselli & Myers 1995), so that hydrodynamic core formation 
models may be more appropriate in this case. Indeed, observations of Class ob- 
jects indicate that protostellar collapse is more dynamic, with supersonic infall 
velocities and large mass accretion rates (> lOCg/G), in cluster-forming clumps 
(e.g. Di Francesco et al. 2001, Belloche et al. 2006). Evidence of coherent, super- 
sonic contraction motions over more than 0.5 pc has even been found in some 
protoclusters (e.g. Motte et al. 2005, Peretto et al. 2006). On small scales, how- 
ever, the compact (~ 0.03 pc) prestellar condensations of the Ophiuchus, Serpens, 
Perseus, and Orion protoclusters are themselves characterized by subsonic lev- 
els of internal turbulence (Myers 2001, Andre et al. 2007) reminiscent of the 
thermal cores of Taurus. Furthermore, the condensation-to-condensation velocity 
dispersion measured in these cluster-forming regions is small and only subvirial 
(Andre et al. 2007). This is consistent with the view that protoclusters often start 
their evolution from "cold" , out-of-equilibrium initial conditions (cf. Adams et al. 
2006), perhaps as a result of external perturbations (cf. Nutter et al. 2006). The 
observed small relative velocity dispersion also implies that collisions between 
condensations in a low-mass protocluster such as LI 688 are relatively rare and 
that, in general, prestellar condensations do not have time to interact with one 
another before evolving into pre- main sequence objects (Andre et al. 2007). Fur- 
thermore, in such protoclusters, the mass accretion rate expected from competi- 
tive, Bondi-like accretion of background gas onto the condensations (e.g. Bonnell 
et al. 2001) is estimated to be at least a factor ~ 3 lower than the mass infall rate 
resulting from gravitational collapse at the Class and Class I stages (Andre et 
al. 2007 - see also Krumholz et al. 2005). Therefore, competitive accretion cannot 
play a dominant role once individual protostellar collapse sets in. On the other 
hand, Bondi-like accretion of unbound gas is more effective before protostellar 
collapse, and may possibly govern the growth of starless, self-gravitating conden- 
sations initially produced by gravitational fragmentation (cf. § 13.11 and Fig. 
toward a Salpeter-like IMF mass spectrum (cf. Myers 2000, Basu & Jones 2004, 
Clark & Bonnell 2005). 

4. Collapse and subfragmentation of prestellar cores 

4.1. Core collapse models: Thermodynamics 

The evolution of gravitationally collapsing cores and the formation of protostars 
are radiation-magnetohydrodynamical processes. Although we need to model 
these processes by solving equations of radiative transfer and magnetohydro dy- 
namics simultaneously in multi-dimensions, a direct calculation of all the equa- 
tions remains challenging. The most sophisticated theoretical models so far are 
either non-magnetic radiation hydro dynamical calculations based on the (flux- 
limited) diffusion approximation or magnetohydro dynamical calculations with 
some prescribed equations of states. 

Here, we first explain the thermodynamics of gravitational collapse revealed by 
dynamical modelling with detailed radiative transfer in spherical symmetry (Lar- 
son 1969; Narita, Nakano, & Hayashi 1970; Winkler & Newman 1980; Stabler, 
Shu, & Taam 1980; Masunaga, Miyama, & Inutsuka 1998; Masunaga &: Inutsuka 
2000a). The classical results based on gray-approximation were confirmed by re- 
cent work that solved the frequency-dependent RHD equations by the "Variable 
Eddington Factor Method" (Masunaga & Inutsuka 2000a). The latter provides 
the time evolution of the apparent spectrum of the radiation field (Spectral En- 
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Figure 9. Temperature evolution at the center of a gravitationally collapsing cloud 
obtained by Masunaga & Inutsuka (2000a) in their radiation hydrodynamical calculation 
of protostellar collapse in spherical symmetry. The first collapse phase corresponds to the 
formation of the first protostellar core that consists mainly of hydrogen molecules. The 
dissociation of hydrogen molecule triggers the second collapse that eventually produces 
the second core^ i.e., a protostellar object. Each of these phases in the temperature 
evolution is characterized by a distinct value of the effective ratio of specific heats, 7cfi. 

ergy Distribution - SED) in addition to the detailed dynamical evolution of the 
protostar. 

Figure [9] shows the time evolution of the central temperature (as a function of 
the central density) in a collapsing cloud. Once compressional heating dominates 
over radiative cooling, the central temperature increases gradually above the low 
value (~ 10 K) found in molecular clouds. At first, the slope of the tempera- 
ture curve corresponds to a ratio of specific heats 7efr = 5/3: T{p) oc p^^^ for 
lOK < T < 10^ K. This apparently monoatomic gas behavior is due to the fact 
that the rotational degree of freedom of molecular hydrogen is not excited in this 
low temperature regime {E{J = 2 — 0)/A;b = 512 K). When the temperature rises 
above ~ 10^ K, the slope becomes that of diatomic molecules (7efr = 7/5). In 
both cases, the effective ratio of specific heats is larger than the critical value for 
gas pressure support against self-gravity: 7efr > 7crit = 4/3. Thus, the col- 
lapsing velocity is decelerated and forms a shock at the surface of quasi-adiabatic 
hydrostatic object, the "first core". Its radius is about 1 AU in spherically sym- 
metric calculations (but is larger by an order of magnitude in 2D / 3D calculations 
with rotation). It mainly consists of II2. The increase in density and tempera- 
ture inside the first core is slow but monotonic. When the temperature becomes 
> 10^ K, the dissociation of H2 starts. The binding energy of H2 is about 4.5 eV 
which is much larger than the thermal energy per hydrogen molecule in this 
temperature regime. Therefore the dissociation of H2 acts as an efficient coolant 
of the gas, which reduces the effective ratio of specific heats below the critical 
value (7ofr < 4/3), and triggers the second dynamical collaps^. In this "second 
collapse" phase, the collapsing velocity becomes very large and engulfs the first 
core. As a result, the first core lasts only for ~ 10^ yr. In the course of the second 
collapse, the central density reaches a stellar value (/?* ~ 1 gcm~^), and a truly 
hydrostatic protostellar object forms in the center. 

t Note that the role of H2 dissociation in this second collapse is analogous to that of the photo-dis- 
integration of Fe as the cause of the pre-supernova gravitational collapse. 
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Radiation hydrodynamic (RHD) calculations automatically produce the time 
evolution of the accretion luminosity and SED of the collapsing object (cf. Ma- 
sunaga k, Inutsuka 2000a) . The resultant luminosity evolution has a sharp growth 
at the formation of the first core and a peak around the formation of the second 
protostellar core in the case of dynamical initial conditions, while it only shows 
a gradual growth in the case of hydrostatic equilibrium initial conditions. This 
difference in the time evolution of the accretion luminosity may provide a useful 
observational diagnostic. 

Molecular emission line profiles of various important species were also calcu- 
lated in self-consistent dynamical models (Masunaga &: Inutsuka 2000b). Re- 
cent three dimensional modelling of protostellar radiation hydrodynamics can 
be found, e.g., in Whitehouse & Bate (2006), Stamatellos et al. (2007b), and 
Krumholz et al. (2007). 

Further evolution includes the T-Tauri phase on Hayashi tracks, for which the 
relevant (Kelvin-Helmholtz) timescale is about two orders of magnitude larger 
than the dynamical timescale (~ 10^ yr) of protostellar collapse, which is only 
accessible by steady state calculations (e.g. Chabrier & Baraffe 2000). 

4.1.1. Dynamical roles of the first and second protostellar cores 

The formation of the 'short-lived' first core plays an important role in the 
dynamical evolution of protostellar collapse, as well as the second core does. 
This is because the associated pressure support helps the formation of a rapidly 
rotating disk-like structure around the central object in the presence of non-zero 
initial angular momentum (Saigo, Matsumoto, & Hanawa 2000). 

Rotationally-driven fragmentation into multiple systems: 

The mechanisms responsible for the formation of multiple systems have been 
the subject of extensive theoretical work (see, e.g., Goodwin et al. 2007 for a 
review). The essence of the results may be summarized as follows. First, gravita- 
tional sub-fragmentation within a collapsing cloud core is not expected during the 
early isothermal phase of evolution, unless the mass of the core is much larger than 
the initial Jeans mass. This result has been found in numerical simulations and is 
supported by semi-analytic calculations in the case of an initially uniform cloud 
with solid rotation (Tsuribe k, Inutsuka 1999a,b). Initial central concentration in 
the density profile makes core fragmentation even more difficult. However, if the 
evolution is followed to the higher density regime where the gas becomes adia- 
batic, a disk-like structure forms which allows another mode of binary formation 
to develop, i.e., disk fragmentation around the central protostar. For example, 
calculations based on a piecewise polytropic equation of state show that the cen- 
tral portion of a collapsing core becomes adiabatic and forms a disc-like structure 
around the central object, which subsequently fragments into "satellite" objects 
(Matsumoto & Hanawa 2003) . The study of realistic radiative transfer effects on 
these modes of binary fragmentation remains an important task for future work. 
The above phenomena correspond to the fragmentation of/around the first core, 
which might account for the formation of binaries with separations ^ lAU. Ob- 
viously formation of binaries with even shorter separations is also expected in the 
disk around the second protostellar core. These multiple epochs of core fragmen- 
tation may result in multiple peaks in the separation distribution of binary stars 
as proposed by Machida et al. (2007c). 

The effects of initial turbulence on binary fragmentation have also been studied 
extensively (see, e.g., Goodwin et al. (2007). Purely hydrodynamic SPH simu- 
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lations of rotating cloud core collapse show that a very low level of initial core 
turbulence (e.g. -Eturb /-E'grav ~ 5%) leads to the formation of a multiple system 
(Goodwin et al. 2004; Hennebelle et al. 2004). In such hydrodynamic simulations, 
fragmentation is driven by a combination of rotation/turbulence and occurs in 
large ( ^ 100 AU) disk-like structures or "circumstellar accretion regions" (CARs 
- cf. Goodwin et al. 2007). These CARs are highly susceptible to spiral instabil- 
ities which always fragment them into small- multiple systems with N > 2 and 
typically N ~ 3-4 within a radius ~ 150 AU (Goodwin et al. 2004). However, 
recent MHD simulations of magnetized core collapse (Price &: Bate 2007, Hen- 
nebelle & Teyssier 2008) show that the presence of an even moderate magnetic 
field strongly modifies angular momentum transport during collapse and at least 
partly suppresses core fragmentation. Therefore, it is unclear at the present time 
whether the collapse of an individual prestellar core typically produces one, two, 
or more stars. 

Implications for the CMF-IMF connection: 

If each prestellar core is the progenitor of A^ ~ 2 — 4 stars, then the IMF of 
individual stars cannot be the direct product of the CMF but results instead 
from the convolution of the CMF with the typical distribution of object masses 
produced by binary fragmentation for one mass of core (cf. Delgado-Donate et 
al. 2003, Goodwin et al. 2008). For realistic binary fragmentation scenarios, the 
IMF will still follow the CMF at the high-mass end (because the majority of each 
core's mass can still end up in one stellar component), but may differ substantially 
from the CMF at the low-mass end. Such a picture is consistent with present 
determinations of the CMF (see § [2]). However, it is presently unclear whether 
the collapse of an individual prestellar core typically produces one, two, or more 
stars. Accordingly, the origin of the low- mass end ( ^ 0.1 Mq) of the IMF is 
highly uncertain. Observationally, this is an area where the future large millimeter 
interferometer ALMA will yield key progress. 

Driving MHD outflovi^s: 

Another effect of the formation of the first core occurs in the MHD evolution 
of the self-gravitating collapsing core. A magnetically supercritical core whose 
rotation axis is parallel to the mean direction of the magnetic field lines leads 
to self-similar collapse as long as the equation of state is isothermal (Basu & 
Mouschovias 1994). Once the first core is formed, a rapidly rotating disk- like 
structure develops owing to the change in the effective equation of state, and its 
rapid rotation winds up the field lines creating a significant amount of toroidal 
magnetic field. This enhanced toroidal magnetic field produces a bipolar outflow 
driven by magnetic pressure (Tomisaka 2002). Thus, the formation of the first 
hydrostatic core plays a critical role in launching the protostellar outflow. A sim- 
ilar process happens again around the second core, with higher ejection velocities 
reminiscent of the observed optical jets and high- velocity neutral winds (Machida 
et al. 2006, 2007a, 2007b). 

The observational detection of the first core would not only confirm the predic- 
tions of RHD modelling, but would also set strong constraints on MHD models 
of protostellar outflows as described in the next section. 
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4.2. Resistive MHD effects and onset of outflows 

4.2.1. MHD modelling with resistivity 

The full modelling of magnetohydrodynamical processes in star formation must 
include non-ideal MHD effects (e.g. Nakano et al. 2002, Tassis &: Mouschovias 
2007). Ambipolar diffusion is important at early times during the low-density 
core formation phase, but Ohmic dissipation is more important at later times in 
the high-density collapse phase. The Hall current term can also be important in 
an intermediate regime ( Wardle 2004) . To account for the dissipation of magnetic 
fields during the formation of protostars, Machida et al. (2006, 2007a, 2007b) used 
the resistive MHD equation with prescribed resistivity in their three-dimensional 
nested grid code simulations. Their basic equations are as follows: 

1^ + V • {pv) = 0, (4.1) 

dv 1 
p— + p{v ■ V)v = -VP - —B X (V X S) - pV</), (4.2) 
ot 47r 

dB 

= \/ X (v X B) +r]V^B, (4.3) 

at 

= 4TrGp, (4.4) 

where p, v, P, B, r/, and <j) denote the density, velocity, pressure, magnetic flux 
density, resistivity, and gravitational potential, respectively. Machida et al. esti- 
mated the resistivity -q in equation (j4.3p according to Nakano et al. (2002) and 
assumed that r] was a function of density and temperature. They further as- 
sumed a barotropic equation of state to mimic the temperature evolution shown 
in Fig. [9l Hence r/ could be expressed as a function of density only: i] = Crj ??o(/o), 
where ?/o(p) is a function of the central density. The initial conditions adopted 
by Machida et al. correspond to a spherical cloud with a critical Bonnor-Ebert 
density profile having a central (number) density pc,o = 3.841 x 10~^^ g cm~^ 
{ficfl = lO^cm"^). In this case the critical Bonnor-Ebert sphere radius, Rc = 
6.45 Cs[47rGpBE(0)]"^/^ corresponds to Rc = 4.58 x 10^ AU. The total mass inside 
the critical radius was Mq = 7.6 Mq. Initially, the cloud was in solid body rotation 
around the z-axis (at a rate Qq) and had a uniform magnetic field (-Binit = 17/iG) 
parallel to the z-axis (or rotation axis). To promote contraction, the density was 
increased by 70% starting from the critical Bonnor-Ebert sphere. 

The various models investigated by Machida et al. can be characterized by a 
single non-dimensional parameter uj, related to the cloud's initial rotation rate, 
and defined using the central density po as co = r2o/(47r Gpo)^''^-The parameter 
UJ, the initial magnetic field strength, the initial angular velocity i^o, the ratio of 
the thermal (ao) and rotational (/3o) energies to the gravitational energy, the final 
magnetic field strength at the center, and the rotation period of the protostar at 
the final snapshot of the calculation are summarized in Table [21 where SR, MR, 
and RR stand for (initially) slow, medium, and rapid rotator, respectively. 

4.2.2. Protostellar outflows & jets 

As a result of their calculations, Machida et al. found that two distinct fiows 
(low- and high-velocity fiows) are driven by the first and second cores. They 
proposed that the low-velocity flow from the first core corresponds to observed 
molecular outflows, while the high-velocity flow from the protostar corresponds 
to observed optical jets. As an illustration of their simulations, a snapshot of 
Model RR is shown in Fig. [TUJ 
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Table 2. Model Parameters and Results f 



Model u 




no [s-i 




ao 


do 


Sfin (kG) 


^fin (day) 


SR 0.003 


17 


7.0 X 10- 


16 


0.5 


3 X 10-^ 


2.18 


3.0 


MR 0.03 


17 


7.0 X 10- 


15 


0.5 


3 X 10-3 


0.40 


2.1 


RR 0.3 


17 


7.0 X 10- 


14 


0.5 


3 X 10-1 







f Representing the thermal, rotational, and gravitational energies as U , K, and W, the relative factors 
against the gravitational energy are defined as ao = U/\W\, and I3q = K/\W\. 

The results of Machida et al. show that the flow driven by the first core has a 
slow speed and a wide opening angle, while the flow driven by the protostar has a 
high speed and a well-collimated structure. The flow speed roughly corresponds 
to the escape speed of the driving object. The difference in the depth of the 
gravitational potential between the first and the second core therefore causes the 
difference in flow speed. Typically, observed molecular outflows and optical jets 
have speeds fout,obs — 30kms"-^ and fjet,obs ~ 100kms~^, respectively. At the 
end of the calculations, the low- and high-velocity flows of Machida et al. only 
have speeds wlvf — 3kms~^ and fnvF — SOkms"^, respectively. However, the 
first and second cores only have masses Mfirstcore = 0.01 and Msecondcore — 
1O~^M0, respectively, at the end of the calculations. Each core grows in mass 
by at least 1-2 orders of magnitude in the subsequent gas accretion phase. Since 
the escape speed increases as the square root of the mass of the central object at 
a fixed radius, the speeds of the low- and high-velocity flows may increase by a 




Figure 10. Bird's-eye view of model RR {I = 12). The structure of high-density region 
[n > lO^^cm-^; iso-density surface), and magnetic field lines (black-and-white stream- 
lines) are plotted. The structure of the outflow is shown by the iso- velocity surface inside 
which the gas is outflowing from the center. The density contours and velocity vectors 
(thin arrows) on the mid-plane of x = 0, y = 0, and z =0 are, respectively, projected in 
each wall surface. 
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Figure 11. Schematic picture proposed by Machida et al. (2007b) for the jet and 
outflow driven from the protostar and the fist core, respectively. 

factor of ~ 10, and reach ulvf — 30kms~^ and uhvf = 300kms~^, respectively, 
which correspond to typical observed values. 

The difference in coUimation between the two flows is caused both by different 
configurations of the magnetic field lines around the driving object and differ- 
ent driving mechanisms (see also Banerjee &: Pudritz 2006 and Hennebelle &: 
Fromang 2008). The magnetic field lines around the first core have an hourglass 
configuration because they converge to the cloud center as the cloud collapses, and 
Ohmic dissipation is ineffective before first core formation. Furthermore, the fiow 
emerging near the first core is mainly a disk wind driven by the magnetocentrifu- 
gal wind mechanism. The centrifugal force is dominant in the low-velocity flow, 
whereas near the protostar, the magnetic field lines are straight, and the mag- 
netic pressure gradient mechanism is more effective in driving the high-velocity 
fiow. The magnetic field lines are stretched by the magnetic tension force near the 
protostar because the magnetic field is decoupled from the neutral gas. However, 
the magnetic field lines are strongly twisted in the region in close proximity to the 
protostar, where the magnetic field is coupled with the neutral gas again. Thus, 
the strong toroidal field generated around the protostar can drive a high- velocity 
fiow, which is guided by the straight configuration of the magnetic field. 

Figure [11] summarizes the main features of the outflows and jets modelled by 
Machida et al. (2007b). Note that the calculations of Machida et al. cover only 
the very early phase of protostar evolution. Further longer-term calculations are 
needed to better understand the correspondence between the models and observed 
protostellar outflows. 

4.2.3. Effects of outflows/jets on star-forming cores/clouds 

In principle, outflows and jets from new-born stars may have strong dynamical 
effects on their environments. It has been proposed that the violent impact of 
fast outflows may result in the destruction of the cores from which the driving 
protostars are born (Nakano et al. 1995, Matzner & McKee 2000), or even in the 
dispersal of the parental molecular clouds. In particular, the local star formation 



24 



efficiency at the level of individual cloud cores is most likely controlled by the 
effects of protostellar outflows: The model calculations of Matzner & McKee 
(2000) give ecore = M^^/Mcore ~ 25 — 75% for various degrees of core flattening and 
magnetization. On a more global level, protostellar outflows may also significantly 
reduce the fraction of total molecular cloud mass going into stars. Thus, this 
process provides an interesting possibility for explaining the observed low star 
formation efficiency in the Galaxy. 

Li & Nakamura (2006) proposed that protostellar outflows could also sustain 
a high level of turbulence in cluster-forming clouds, but Banerjee et al. (2007) 
argued that this was actually difficult. Obviously, a quantitative examination of 
all these processes will require appropriate dynamical modelling of protostellar 
outflows, which remains to be done both theoretically and observationally. 

5. Conclusions: Proposed view for the star formation process 

The prestellar core mass function (CMF) appears to be consistent with the 
stellar IMF between ~ 0.1 Mq and ~ 5Mq, although large uncertainties remain 
especially at the low- and high- mass ends (cf. § ED- Small internal and relative 
motions are measured for these prestellar cores, implying that they are much less 
turbulent than their parent cloud and generally do not have time to interact before 
collapsing to (proto)stars (cf. ^ I3.2.5p . These results strongly support scenarios 
according to which the IMF is largely determined at the prestellar stage. 

None of the extreme scenarios proposed for the formation of prestellar cores can 
explain all observations. Pure ambipolar diffusion is too slow to be the main core 
formation mechanism, for typical levels of cloud ionization. Purely hydrodynamic 
pictures have trouble accounting for the inefficiency of core formation and the 
detailed velocity structure of individual cores. A mixed scenario (cf. Nakamura &; 
Li 2005, Basu et al. 2008) may be the solution: supersonic MHD turbulence in a 
molecular cloud close to magnetic criticality generates seed cores, which grow in 
mass until they become gravitationally unstable and collapse in a magnetically- 
controlled fashion while decoupling from their turbulent environment. Such a 
mixed picture has several advantages. On the one hand, relatively strong magnetic 
fields prevent global collapse and lead to inefficient core/star formation on large 
(GMC) scales. On the other hand, the rate of ambipolar diffusion is enhanced 
by turbulence in shock-compressed regions (e.g. Zweibel 2002, Fatuzzo & Adams 
2002, Nakamura & Li 2005). The difference between the clustered and the isolated 
or distributed mode of star formation may result from local variations around 
the critical mass-to-fiux ratio and/or from the presence or absence of external 
triggers. 

Cluster-forming clumps such as L1688 in Ophiuchus or NGC2264-C in Mono- 
ceros are likely (slightly) magnetically supercritical. There is some evidence that 
these clumps are in a state of global collapse induced by large-scale external 
triggers (e.g. Peretto et al. 2006, Nutter et al. 2006). In this case, the local star 
formation efficiency within each condensation is high ( ^ 50% - Motte et al. 
1998) and a promising core formation mechanism is purely gravitational, Jeans- 
like fragmentation of compressed cloud layers (cf. Palous 2007, Peretto et al. 
2007), possibly followed by subsequent core growth (e.g. Basu &: Jones 2004). 

Regions of more distributed and less efficient star formation, such as Taurus 
or the Pipe Nebula, are likely slightly subcritical or nearly critical. The size of 
individual cores is then larger (e.g. Motte et al. 1998), the local star formation 
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efficiency is lower (~ 15%-30% - Onishi et al. 2002, Alves et al. 2007), and the 
feedback from protostellar outflows may be more important in limiting accretion 
and defining stellar masses (e.g. Shu et al. 2004). 

To fully understand how stars form and how the IMF comes about, a better 
knowledge of the presently poorly known initial conditions for molecular cloud 
formation is required, a subject of growing interest (see, e.g., Hennebelle et al. 
in this volume). It is also crucial to further investigate the processes by which 
prcstcUar cores form, evolve, and eventually collapse and fragment into multiple 
protostellar systems. With present submillimctcr instrumentation, observational 
studies are limited by small-number statistics and restricted to the nearest re- 
gions. The advent of major new facilities in the coming years should yield several 
breakthroughs in this field. With an angular resolution at 75-300 fim compara- 
ble to, or better than, the largest ground-based millimeter-wave radiotelescopes, 
Herschel, the Far InfraRed and Submillimeter Telescope to be launched by ESA 
in 2008 (cf. Pilbratt 2005), will make possible complete surveys for prestellar 
cores down to the proto-brown dwarf regime in the cloud complexes of the Gould 
Belt (cf. Andre & Saraceno 2005). High-resolution (0.01" -0.1") studies with the 
'Atacama Large Millimeter Array' (ALMA, becoming partly available in 2011, 
fully operational in 2013 - cf. Bachiller 2008) at ~ 450 jim - 3 mm will allow us to 
probe the kinematics of individual condensations in distant, massive protoclus- 
ters. Complementing each other nicely, Herschel and ALMA will tremendously 
improve our global understanding of the initial stages of star formation in the 
Galaxy. 
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